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t~^ ■ Abstract: We present an implementation of the vector boson pair production processes 

ZZ, W^W^ and W Z within the POWHEG framework, which is a method that ahows the 
interfacing of NLO calculations to shower Monte Carlo programs. The implementation is 
built within the POWHEG BOX package. The Z/"f* interference, as well as singly resonant 

k><( \ contributions, are properly included. We also considered interference terms arising from 

^ \ identical leptons in the final state. As a result, all contributions leading to the desired 

four-lepton system have been included in the calculation, with the sole exception of the 
interference between ZZ and W^W~ in the production of a pair of same-flavour, oppositely 
charged fermions and a pair of neutrinos, which we show to be fully negligible. Anomalous 
trilinear couplings can be also set in the program, and we give some examples of their effect 
at the LHC. We have made the relevant code available at the POWHEG BOX web site. 
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1. Introduction 

The pair production of electroweak vector bosons is of great interest to the particle physics 
community. The leptonic decay of electroweak boson pairs has been intensively studied at 
the Tevatron [1-9], while first measurements of W^W~ production were recently published 
by both ATLAS [10] and CMS [11]. 

Vector boson pair production is interesting for various reasons. First of all, the process 
is interesting in itself, as a test of the non-Abelian nature of the electroweak force. Any 
deviation from the Standard Model tri-vector boson couplings would indicate the presence 
of new physics. Electroweak boson pair production is also an irreducible background to 
moderately heavy Higgs production, where the Higgs decays into electroweak bosons. It is 
therefore essential that we are able to make accurate predictions for these processes. 

In general, calculations at leading-order (LO) in perturbative QCD (pQCD) have a 
large dependence on the unphysical factorization and renormalization scales. In order 
to limit these uncertainties, it is necessary to extend the pQCD calculations to next- 
to-leading order (NLO). Analytic formulae for the LO and NLO matrix elements for a 
weak boson pair decaying to leptons, written succinctly in the helicity formalism, were 
presented in ref. [12] ^. These have been implemented in publicly available programs, such 
as MCFM [14,15], enabling fast computations for these processes. 

However, the presence of soft and collinear divergences in the final state means that 
any computation performed at NLO can only provide accurate predictions for inclusive 



^Although qq — )> W^W with decays was previously considered to 0{as) in Ref. [13], the virtual 
correction to the spin correlations were not computed there. 



quantities. On the other hand, exclusive quantities are often of great interest when an- 
alyzing an event - for example, any quantity including a jet variable in weak boson pair 
production. In order to compute exclusive quantities accurately, parton shower programs 
must be used. Two methods exist to interface NLO results to parton shower programs: 
MCQNLD [16] and POWHEG [17,18]. In fact, W^W^ -pair production was the first process 
studied in the former [16], while ZZ production was the first practical demonstration of 
the latter approach [19]. Furthermore, W^W^ , W Z and ZZ production have been 
studied within the POWHEG framework in refs. [20,21]. Notwithstanding this, a full, pub- 
lic NLO implementation of electroweak boson pair production, including leptonic decays, 
Z/7* interference and non-resonant graphs, is still missing, and is badly needed by the 
experimental collaborations in order to reliably simulate their backgrounds away from the 
resonant region. The purpose of the present work is to fill this gap by providing NLO 
implementations for weak boson pair production that include all diagrams that lead to 
the desired 4-lepton final states. The only effect not included is the interference between 
the W~^W^ and ZZ pair production when the final state consists of two opposite charged 
leptons and a neutrino-antineutrino pair of the same family. We show, however, that this 
effect is fully negligible. 

A formally NNLO contribution to these processes, initiated by gluon fusion, is also 
known to be important. We do not include it in our program. Calculations of these 
processes are already available in the literature [15,22-24], and in particular the gg2ZZ 
and gg2WW generators of refs. [22-24] can be easily interfaced to shower program, and 
are currently used by the experimental collaborations for this purpose. 

The rest of this paper is organized as follows. In section 2, we describe the NLO 
calculation of the vector boson pair production processes, and the implementation of these 
processes in the POWHEG BOX framework [25]. We present our results in section 3. The effects 
of the POWHEG implementation and subsequent parton showering, common to all three vector 
boson pair production processes, are demonstrated for ZZ production. We also discuss the 
effects of anomalous trilinear boson couplings in W^W~ and WZ production, focusing on 
the potential of the LHC to improve on the existing bounds on these couplings. In section 
4 we give our conclusions. In the two Appendices we discuss some technical details. 

2. Description of the implementation 

We took the matrix elements for the Born, real and virtual amplitudes from the calculations 
of Ref. [12], borrowing heavily from their implementation in the MCFM package [14, 15]. At 
variance with the MCFM package, however, we have used the Breit-Wigner form for the 
off-shell vector bosons propagators 

z (2 I) 

Sy-M^+iFyMy' ^ ' ' 

rather than the Baur-Zeppenfeld form [26]. We can also enforce the use of the running 
width variant (see manuals for details) 

1 



Sy-M^+ iVySy/My 



2.2 



> 

CD 
O 


14 
12 


£ 


10 



+ 



8 
6 


T3 


4 




2 


^ 


1.1 


O 


1.0 



Interference 

No Interference 

WW 

ZZ 




i^jfe ii ^jfa ^ g i 3ii ^ \j j Efe^6j„t^ i j^j^ i ^ff^jai} i ^ i ^ i j^^jJ i j'Wj^^ 



1 0.9 

20 40 60 80 100 120 140 160 180 200 

Figure 1: Distribution of Afe+g- in the process pp —> e^e^VeVe at leading order and centre of mass 
energy of 7 TeV. This demonstrates the negligible effect of interference between W~^W~ and ZZ 
production of this final state. The solid green and magenta lines show the separate contribution 
from W^W" and ZZ production respectively. 



We have seen, however, no appreciable differences when this last form is used. 

Z/7* interference effects are included. Singly resonant amplitudes, which are essential 
if we wish to describe kinematic regions where only one of the vector bosons is resonant [15] , 
are also included. Furthermore, we include interference terms when there are identical lep- 
tons in the final state, as in ZZ — t- e^e^e^e^, and ZW^ — t- e~e^e^v. These effects are 
small, but become relevant when all oppositely charged lepton pairs are far off-resonance.^ 
Interference terms of this sort never arise in VF"^W^~-production. However, W^W~ and 
ZZ production may lead to identical final states, such as e~e+z^eJ^e- We have studied the 
corresponding interference effect at leading order with Madgraph [27] . As shown for exam- 
ple in Fig. 1 in the case of the Mg+g- distribution, such effects are completely negligible, 
and are therefore not included in the following. Thus, if one wishes to investigate leptonic 
final states like e~e^V(,i>^^ both a W^W" and a ZZ sample should be generated separately. 

The construction of the Born phase space requires some care, in order to efficiently 
probe the resonance regions. In Appendix A we discuss this problem with some detail. Fur- 
thermore, as in the case of W production, a problem due to vanishing Born configurations 
has to be dealt with. This is explained in Appendix B. 

The three programs perform quite differently. The W Z program can complete the 
preparation stage in roughly one hour, involving about 6 million calls, and can generate 
one million events (at the Les Houches level) in six hours. The WW program has similar 
performance, but the preparation stage may require twice as much. The ZZ program is 
much slower. Even without including interference for identical fermions, the preparation 
stage takes of the order of 6 hours, and one million events are generated in 40 hours. It 



^Since interference effects are very small, and their calculation more than doubles the computing time, it 
is possible to switch them off by setting an appropriate flag in the powheg. input file. Similarly, it is possible to 
generate on-shell vector bosons in the zero-width approximation or to remove the single resonant diagrams, 
as illustrated in the implementation manuals. 



is thus convenient, in this case, to use the POWHEG BOX facihties to run the program in 
parallel. Requiring interference costs more than a factor of two in execution time. The 
poorer perfomance of the ZZ program is easily tracked back the the presence of a larger 
number of helicity configurations that contribute to the process. 

Several checks have been performed to test the correctness of the implementation. 
The Born and real matrix elements have been checked against MadGraph [27] for an arbi- 
trary phase-space point. The ratio of the residues of the single and double poles of the 
virtual amplitude to the Born amplitude were compared against the analytically known 
expressions. The finite part of the virtual amplitude was checked against an independent 
Feynman diagram based program, developed by some of the authors, which uses the OPP 
subtraction method [28]. The MCFM program implements the same processes, but does not 
include the interference terms for identical fermions. We have thus verified that, when 
the interference terms are not included or not present, the NLO output of the POWHEG im- 
plementation, which can be obtained by setting an appropriate option in the program, is 
consistent within errors to that of MCFM for a large number of distributions. Furthermore, 
the POWHEG output at the level of the Les Houches Interface events has been compared for 
on shell Z bosons to the previous implementation of ref. [19]. Full agreement was found, 
thus demonstrating that the mechanism used by POWHEG to generate radiation operates in 
the same way as in the implementation of ref. [19], which in turn was thoroughly compared 
to the MCONLO implementation. 

3. Results 

We now present some results obtained with our POWHEG implementation. Our basic setup 
for the electroweak input parameters is the following. We use M^, M^ and G^ as basic 
parameters of the Standard Model (SM), and define all remaining quantities according to 
the leading order SM relations. Thus we have 

cos6'^ = — — , acm = sm fc*^ (3.1) 

Mz TT 

This corresponds to the so called "G^ scheme", advocated in ref. [29]. We adopt the 
following PDG [30] values for the independent parameters 

Mz = 91.1876 GeV, M^^, = 80.399 GeV, G^ = 1.166364 x 10"^ GeV"^ . (3.2) 

For the width of the W and Z bosons, we take their LO value computed in the above 
scheme, except that we correct the hadronic width with a factor (1 -|- as(Mn/ )/7r) and 
{1 + as{Mz)/7r), respectively. In this way, the branching fractions of the W/Z consistently 
add up to one. For the W and Z width we get 2.0997 GeV and 2.5096 GeV, respectively. 
These values are around half a percent away from the current PDG values [30], 2.085 it 
0.042 GeV and 2.4952 it 0.0023 GeV, respectively, so that, in fact, the measured values 
could be used instead. Details can be found in the file smcouplings.f , which can easily 
be modified by users preferring other ways to implement SM couplings. As default, we set 
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Figure 2: The effect of different levels of approximation of the calculation oi pp — > e+e e+e at 
7 TeV for the invariant mass of the four-lepton system. 



the factorization and renormalization scales equal to the mass of the four-lepton system, 
Mil, and our default PDF set is the MSTW 2008 NLO [31]. 

The three processes we examine are very similar to one another. For each one of 
them we have analyzed several distributions, such as: the rapidity, pseudorapidity, trans- 
verse momentum and invariant mass of each lepton, each combination of lepton-lepton and 
lepton-antilepton pairs, and the four lepton system; the distance in rapidity, pseudorapid- 
ity, azimuth and AR distance (where AR = y A(^2 + 6r]'^) of all lepton-antilepton pairs; 
and the rapidity distribution of the hardest jet and its distance in rapidity from the four 
lepton system, with several transverse momentum cuts. For each distribution, the pure 
NLO result has been compared to the POWHEG result before the shower. All distributions 
behave as expected. In the present paper we illustrate a selected set of distributions with 
realistic cuts in the ZZ case only. We compare the various levels of approximation that 
one can include in the NLO calculation, and the full NLO result with the event output by 
POWHEG before and after the shower. 

We remind the reader that the POWHEG output can be easily interfaced to all gen- 
eral parton shower event generators [32-35] that comply with the Les Houches event files 
format [36]. In the present work we use PYTHIA version 6.4.25 [32]. 

We illustrate our implementation of W'^W' and WZ production by studying the effect 
that anomalous couplings have on these processes at the LHC. 

3.1 ZZ Production 

We first illustrate the differences that arise at different levels of approximation in our 
generator in the ZZ case. In Fig. 2 we show three predictions for the invariant mass of 
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Table 1: Total cross section for pp -^ ZZ -^ e+e^/i+/i^, with the only cut Mi+i- > 20 GeV 
on both lepton pairs. The central values are for the scales fip = fiR = Mn, and the integration 
error in the final digit is shown in parentheses. The theoretical uncertainty is obtained by doubling 
and halving independently the scales iip and fip, with respect to their central value, excluding the 
combinations that yield ^ip/ Hr — 4 and UnJ Hf = 4. 



the four-lepton system in ZZ production, when both Z bosons decay into electrons. The 
cuts are as documented further on, but for the purpose of understanding the present plot, 
the only relevant cut is the 20 GeV one on the invariant mass of opposite sign leptons. 
The curve labelled as (c) is the full calculation of the four lepton production process, and 
thus is the more accurate result that our generator can provide. In (b), we suppress the 
interference effects due to identical fermion pairs in the final state. The ratio of (c) over 
(b) is presented in the lower panel. We see that there are sensible differences around the 
Z peak. This is easily understood. The singly-resonant production mechanism includes 
the production of a single Z boson that decays into four leptons, giving rise to the Z peak 
observed in the figure. In this decay process both lepton pairs are off resonance, and thus 
we can expect non-negligible interference effects. The knee visible in the figure at roughly 
M^i « 110-120 GeV is due to the onset of the production of an on-shell Z decaying into 
a lepton pair, with the second lepton pair arising from a Z* /^* decay (this mechanism is 
displaced to 110-120 GeV because of the 20 GeV cut on the minimum invariant mass for 
a lepton pair). In this kinematic regime, with an on-shell Z decaying into a lepton pair, 
the interference is suppressed by the phase space. Finally, the curve labelled (a) does not 
include single resonant graphs or interference for identical fermions. It is in good agreement 
with the full result only if we are sufficiently beyond the threshold for the production of one 
on-shell Z boson. The approximation (b) is presently implemented in the MCFM package. 
The approximation (a) is implemented, for example, in PYTHIA, where off-shell effects and 
Z/^* interference are accounted for, but singly resonant graphs are not included. Because 
of the small differences between the single resonant (b) and the full calculation (c), for 
simplicity we proceed by considering decays to e^e^/U+^u^ only. 

In the following, we employ a cut of Mi+i- > 20 GeV for each electrically charged, 
same-flavour lepton-antilepton pair l^l^ , which ensures that divergences due to low-virtuality 
photon emissions are avoided. We remark here that, if interference effects due to identical 
fermions are neglected, the same flavour and different flavour cross sections can be simply 
related; for example, the e~^e~ fi^fi^ cross section is twice the e^e~e'^e~ one. We should 
however remember that cuts applied to the final state are generally different: in the case 
of e~^e~e~^e~ we limit the invariant mass of all four possible choices of oppositely charged 
lepton pairs, whereas in the e~^e~n^n~ case there are only two such pairs. Thus, the cross 
sections are related by a factor of 2 (if we neglect interference) only if we impose additional 



mass cuts on the e^/x^ and e^/U+ pairs in the e^e^ fi^n^ case. 

In Table 1, we show the ZZ — t- e~e^ fi~ h~^ cross-section, at leading and next-to-leading 
order, and show the effects of varying independently the renormalization and factorization 
scale by a factor of two in either direction, excluding the values that lead to a ratio of the 
two scales equal to 4 or 1/4. We thus consider the seven combinations: fip = fi^ = M41, 
2iip = fijf = M41, p.p/2 = Hr = M41, lip = 2fj,ji = Mil, fJ-F = I^r/2 = M41, 2fj,p = 2^^ = 
M41, fip/2 = ijLr/2 = M41, with the first one taken as the central value. This is done for 
three PDF sets: MSTW08 [31], CTIO [37] and NNPDF2.1 [38]. Note that we employ NLO 
parton distribution functions also for LO cross-sections. The scale uncertainty is just 1-2% 
at leading order, and slightly smaller at next-to-leading order. We observe, however, that 
the NLO result does not lie within the LO scale variation band. The anomalously small 
scale dependence of the LO cross-section is due to the fact that the Born level process is 
of order zero in the strong coupling constant, and thus the only scale dependence comes 
from the PDFs. Furthermore, in pp collisions new quark-gluon channels open up in the 
NLO cross-section that compete sensibly with the LO one. Notice also that the sfg-initiated 
contribution, which is formally of NNLO (and thus is not included in the table), is larger 
than the scale variation seen here. All these considerations lead to the conclusion that 
scale variation alone underestimates higher order effects for this process. The difference in 
cross-sections between different PDF sets is similar in size to the scale variation. 

Experience with the POWHEG implementation of other processes has taught us that, for 
processes that involve only initial state radiation (ISR) , one expects the effect of showers 
to be modest, a fact that we will verify also in the present context. 

We begin by carrying out a detailed comparison of the POWHEG results at the level of 
the Les Houches events output (the LHE level from now on) with the NLO results. We 
impose the following cuts: the highest pt lepton must have pt > 20 GeV, and the other 
three hardest leptons must have pt > 10 GeV. We also require for the lepton rapidities 
that \rji\ < 2.5. Jets are reconstructed using the anti-kt algorithm [39] as implemented in 
Fast Jet [40], with R = 0.6. In Fig. 3, we show results obtained by analyzing the LHE 
files directly, without further showering, compared to the NLO result. Six distributions 
are presented: M^^, M^j, M41, pt,zi, Ve^i ^^^ Uzi, where y is rapidity and M denotes an 
invariant mass distribution. For every phase space point, we assign Zi to the lepton pair 
of the same flavour whose invariant mass is closest to Mz and label the other lepton pair 
as Z2- This explains the difference between the first two plots. In fact, the lepton pair 
closest to the Z mass tends to have a much reduced contribution from 7* production. In 
general, we see good agreement, at the level of few percent, between the two predictions 
for these inclusive quantities. Other inclusive variables exhibit a similar behaviour. 

We now show the transverse momentum distribution of the four lepton system, pt^4i = 
I ^iPt,i\, in the left plot of Fig. 4, which at NLO is identical to the transverse momentum 
distribution of the jet. We see large differences between NLO and LHE results. At low pt^n, 
the LHE result is much smaller than the NLO one. This fact is typically understood as 
originating from the Sudakov suppression which is present in the POWHEG output, but not in 
the NLO result, which diverges at small Pt,4i- At large pt,4i the LHE result overshoots the 
NLO one by about 50%. An enhancement of the large pt^4i tail of the POWHEG distribution 
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Figure 3: Comparison of NLO to LHE events for pp — > e^e^fi^ii^ at 7 TeV for a number of 
observables, for our default set of cuts. The definitions of Zi and Z2 are given in the text. 



is also observed in processes like single vector boson production and gluon-fusion Higgs 
production. The origin of this effect is discussed in [41,42]. In essence, in POWHEG, the large 
transverse momentum distribution in the real emission process is enhanced by the NLO 
K-factor. We thus expect an enhancement of the order of ~1.25 in our case. Ho'wever, 
we see that this effect accounts only for a fraction of the enhancement observed in the 
plot. The remaining enhancement is easily traced to arise from the different choice of 
scales used in the NLO calculation and in the generation of LHE events, the former being 
equal to the invariant mass of the four lepton system, and the latter being instead taken 
equal to its transverse momentum. These points are illustrated in the right plot of Fig. 
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Figure 4: In the left plot, wc compare the transverse momentum distribution of the four-lcpton 
system in the process pp -^ e^e^ fi'^ii^ at 7 TeV, obtained with the fixed order NLO result, and 
at the Les Houches events level. On the right plot the same comparison is carried out for the 
fixed order NLO result (NLO), and the fixed order NLO calculation with the renormalization and 
factorization scale in the real contribution set to pt^Ai (NLO-b), versus a POWHEG result at the Les 
Houches event level with the bornonly option set (LHE-c). In the lower panel, the curve labelled 
b is the ratio of LHE-c/NLO-b. 



4 where we plot the same quantity using the bornonly option in POWHEG^ compared to 
the standard NLO result, and to the NLO result obtained by setting the factorization and 
renormalization scales equal to pt^n in the real contribution. The ratio of the POWHEG result 
with the bornonly option set, over the NLO result with the pt^/^i scale choice, eliminates two 
causes of the difference (i.e. the different scale choice, and the B/B factor), leaving only 
the Sudakov effect, represented in the figure by the ratio of the LHE-c versus the NLO-b 
result (curve marked as b in the lower panel). Observe that the ratio LHE-c/NLO is about 
25% lower than the LHE/NLO ratio, thus showing that the K-factor effect is indeed of 
the order of 25%. Finally, we remark that the scale choice ^ = pt^n is more appropriate 
for the description of the radiation tranverse momentum than our default one in the NLO 
calculation. Thus, in this respect POWHEG does the right thing. 

We now examine, for the same observables shown in Fig. 3, the effect of a subsequent 
shower, obtained by interfacing the POWHEG output with PYTHIA. When showering, we turn 
off electromagnetic radiation.^ The results can be seen in Fig. 5. As expected, we see very 
little effect from the parton shower for these observables. We remark that in POWHEG, by 
construction, the invariant mass of the two- and four-lepton systems are in fact the same 
at the NLO-level, LHE- level, and after parton shower (regardless of the shower program). 

Turning now to exclusive variables, we show in Fig. 6 the distribution of the jet pseudo- 
rapidity r/j, and the rapidity difference between the jet and the ZZ system, yjzz-, comparing 
LHE with both NLO and the effects of showering with PYTHIA. In addition to the lepton 
cuts already described we have jet cuts of ptj > 20 GeV, and |r/j| < 3.5. We see a noticeable 



^This option eliminates the K-factor enhancement by replacing B -> _B in the POWHEG code, see [41,42] 
for details. 

*This is adequate in this work for purpose of illustration, since an experiment will generally correct the 



lepton energy for electromagnetic radiation emitted coUinearly. 
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Figure 5: Comparison of LHE to PYTHIA showered events for pp — > e+e^/i+/x^ at 7 TeV for a 
number of observables, for our default set of cuts. Tlie definitions of Zi and Z2 are given in the 
text. 

difference between NLO and LHE results in this distribution. This is a consequence of the 
jet Pi-cut and of the difference in the jet pt distribution at NLO and LHE level (see Fig. 4). 
However, a subsequent parton shower has a small effect on the LHE results. 

3.2 W+W/WZ Production 

In this Section we present predictions for W~^W~ /W Z production decaying to four leptons. 
As for ZZ production, the NLO calculation includes Z/^* interference, single resonant 
diagrams, and when relevant, interference due to identical leptons. 

In Tables 2 and 3, we show the cross-sections for WZ — t- e^e~/x^i/u and W~^W~ — >• 



10 




Figure 6: The pseudorapidity of the jet rjj (top row) and the rapidity difference between the 
jet and the ZZ system, yjzz (bottom row), for pp — >■ e+e~/^+/x~ at 7 TeV. Additional cuts of 
Ptj > 20 GeV and \qj\ < 3.5 are apphed. The first column shows a comparison between NLO and 
LHE, and the second column a comparison between the pure LHE and the events showered with 
PYTHIA. 
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Table 2: As in table 1 for pp -> W^Z - 
lepton pair originating from the Z boson. 



1.1 v^e'^e , with the only cut Mg+g- > 20 GeV on the 
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Table 3: As in tabic f for pp -^ W+W' 



h'e^J■ t'/j, with no cuts applied. 



/i-i/^e+Pe, for MSTW2008, CTIO and NN2.1 NLO PDF sets. The theoretical error due 
to scale variation is computed with the same method of table 1. For WZ and W^W~ 
production, we find scale uncertainties of around 1% at leading order, and 2-3% at next- 
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to-leading order. Again the NLO result is not contained in the LO scale variation band. 
This pattern was already observed and discussed in the ZZ production case. 

Our implementation also allows one to study effects due to non Standard-Model like 
(anomalous) WWZ and WWj couplings. While the LHC will probe new physics at the 
TeV scale directly, anomalous trilinear gauge couplings (ATGCs) indirectly probe physics at 
larger scales, since they arise when high-energy degrees of freedom are integrated out. Both 
the Tevatron [7,8,43] and LEP [44] were able to place quite stringent bounds on anomalous 
trilinear couplings. However, since their effects are enhanced at high energies, one may 
expect even better bounds from the LHC. Indeed, CMS already presented bounds on the 
anomalous couplings appearing in an effective Lagrangian with the parametrization of 
ref. [45] (HISZ from now on) without form factors [11]. It is then interesting to understand 
when the LHC running at 7 TeV will be in a position to improve on existing Tevatron and 
LEP bounds [46]. 

Following ref. [45,47,48], we parametrize the most general terms for the WWV vertex 
(y = 7, Z) in a Lagrangian that conserves C and P as 



^yyy y yy /ii/yy y j ^ n yy ^yy yv -r ^ " fiu" p 

(3.3) 
where W^,y = d^Wy — d,yW^, Qwwz = — ecot^^y and Qww-^ = — e. In the SM gX = n\ = 1, 
and A^ = 0. Any departure from these values would be a sign of new physics. Therefore 
anomalous couplings are usually written in terms of the deviation from their SM value, e.g. 
^gX = gX ~ ^ ^tc. In our POWHEG generator, all six parameters can be set independently. 
If one imposes gauge invariance under abelian (electromagnetic) gauge transformations, 
the parameter A.gl vanishes. This still leaves five independent anomalous couplings. The 
LEP groups use a parametrization in which the number of independent couplings reduces 
to three. This is a consequence of only including operators of up to dimension six in their 
effective Lagrangian. One can write the couplings in terms of the three parameters Uw, Ciw4> 
and aBih 



Agf = — ^—— , A^ = A^ = a^, AtC = a^d) + o.b6 , Ak^ = a^y^ — tan^ OwOig^ 



(3.4) 



cos^ 6^^ 
Note that this implies the relation 

Ak^ = Afff - AkT tan^ Ow (3.5) 



In the HISZ model [45], it was also suggested that one may set aw(j> = oiscj) as a further 
simplification, leaving only two independent parameters. This modified setup is sometimes 
used in experimental searches. 

In the presence of anomalous couplings, the effective Lagrangian of eq. (3.3) gives rise 
to interactions that violate unitarity at high energy. Thus, in order to achieve a more 
realistic parametrization, the couplings are multiplied by form factors, that embody the 
effects arising from integrating out the new physics degrees of freedom. The precise details 
of the form factors therefore depend on the particular model considered. Paralleling here 
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Figure 7: The pt distribution for the leading lepton in T4^+iy^-production is shown on the left, 
and the integrated cross-section as a function of the cut on the transverse momentum of the leading 
lepton is shown on the right. The results are shown using Standard Model couplings (in red, labeled 
'SM'), as well as with the anomalous couplings of eq. (3.7) (in green, labeled 'ATGC, no FF'). The 
effect of a form factor given in eq. (3.6) is also shown, both for a value of A = 5 TeV and A = 2 TeV 
(in pink and blue respectively). 



the discussion of ref. [49], we assume that all anomalous coupling Ag are modified as 



A5 



A<7 



:i + M2JA2)2 



(3.6) 



where MyY is the invariant mass of the vector boson pair and A is the scale of the new 
physics. 

In the LEP study of ref. [50] the value of one or two of the parameter (s) is fixed to their 
SM values, and determines the range of values of the remaining parameter(s) consistent 
with their data. We choose the maximum deviation of all parameters from their SM values 
allowed by these fits. While it is true that this point in parameter space would lie outside 
the LEP bounds where all three parameters are allowed to vary, we choose it for illustrative 
purposes. 

The values of the parameters are 



A5^ 



-0.027, 



AA^ = AA^ 



-0.044, A/- 



-0.112, 



(3.7) 



and where (using the LEP groups parametrization) Ak^ is obtained from the above using 
eq. (3.5). We note that these values for the ATGCs also fall within the unitarity bounds 
of ref. [51] for A < 7 TeV. In the left-hand plot of Fig. 7 we show the pt distribution for 
the hardest lepton (which is known to be very sensitive to ATGCS) in I^"'"VF~-production 
using the above ATGCs, with different values of the form factor parameter (A = 2 TeV, 
5 TeV and oo). The effect of the anomalous couplings is evident in the high-pt tail of the 
distribution. The form factor damps the effect of the ATGCs. Its impact increases with 
decreasing A, while A = cxd corresponds to no form factor. In the right-hand plot of Fig. 7, 
we show the integrated cross-section as a function of the cut on the pt of the leading lepton. 
Fig. 7 gives an indication of the kind of signal one would see at the LHC. Depending upon 
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Figure 8: As in Fig. 7 ioi W Z production. 



the precise form of the new physics, the signal would lie somewhere between the A = 2 TeV 
and A = oo predictions. We note that in the region pt > pt,min, with 150 < pt^min ^ 250 
GeV, the cross section is strongly affected by the anomalous couplings. Furthermore, its 
value is of the order of 1 to 10 fb, and thus such a difference should be observable with a 
few inverse femtobarns of data. It may then be possible for the LHC to improve on the 
LEP bounds on anomalous trilinear boson couplings already by the end of the year. 

Similar results are obtained in the case of WZ production, and are shown in Fig. 8 for 
W~ Z. This process is sensitive only to the WW Z tri-linear coupling and not to the WW'j 
couplings. In this case the cross section is smaller and one expects to make less stringent 
bounds than in the W^W~ study. 



4. Conclusion 

We have implemented the production of vector boson pairs W~^W~ , W Z and ZZ in the 
POWHEG BOX. Off-shell effects, Z/'j* interference, single resonant amplitudes and interfer- 
ence between same-flavour leptons have all been included. The only effect not included 
is the interference between W~^W~ and ZZ processes, when both decay to same-flavour 
leptons and neutrinos. We have, however, shown that this interference is entirely negligible, 
so that it suffices to consider these processes separately. We have shown that, for variables 
inclusive in the hardest jet (or equivalently in the radiated parton), there is very little 
difference between NLO and LHE results, whereas noticeable differences arise for variables 
that probe the hardest jet. Neither kind of observable is much affected by the subsequent 
parton shower. 

Our POWHEG implementation allows the study of the effect of anomalous trilinear boson 
couplings in WZ and W^W^ production. We demonstrated this for both processes at the 
LHC with y/s = 7 TeV, and looked at the distribution of the pt of the hardest lepton. 
We saw that a set of anomalous couplings allowed by the LEP bounds leads to a large 
deviation from the Standard Model result in the high pt bins of this distribution, to such 
an extent that the LHC may be able to probe this by the end of the present year. 
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We have not included in this work the W'j or Z7 processes with an on-shell photon. 
This is because NLO calculations with final state on shell photons have to dealt with in 
a different way. Photons can be radiated in the collinear direction by light quarks, and 
thus generate singularities that require special attention. Including photons in a PDWHEG 
generator is however possible, and in ref. [52], for example, the 77 production process was 
implemented. We thus defer this problem to future work. 

Finally, the code for our generators has been made available at the POWHEG BOX website 
http://powhegbox.mib. infn.it. 
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A. Phase space sampling 

Although the POWHEG implementation of the processes we consider is relatively straightfor- 
ward, some care is needed in the generation of the Born phase space, to ensure that the 
resonant regions are probed efficiently. Thus, when implementing processes with no lepton 
interference, like ZZ — t- e^e^ fi^fi^ , the phase space generator should include the virtuality 
of the e"'"e~ and /U+^u^ pairs as integration variables, so that the importance sampling can 
be performed near the Z peak or for low invariant mass, where the pole of the photon 
propagator starts to count. In ZZ and WZ production, a cut on the lepton pair mass for 
oppositely charged same flavour leptons must be imposed, in order to stay away from the 
photon pole. 

When lepton interference is present, like in ZZ — )■ e^e~e'^e~, there are two ways to 
assign the leptons to the vector resonance, and using all the corresponding mass combina- 
tions as integration variables becomes too cumbersome. In this case, in order to maintain 
an efficient importance sampling, we do the following. Calling 1 and 2 the two possible 
resonance assignments, we multiply the cross section by a factor 2F2/{Fi + F2), where the 
Fi/2 are factors that suppress the vector poles for the 1 or 2 resonance assignment. This 
is possible since the two resonance assignments differ only by a permutation of identical 
fermions, and the cross section is symmetric under 10 2 exchange. After multiplying the 
cross section by this factor, the poles of the second region are suppressed, and one can 
adopt a phase space importance sampling as if there was only the resonance assignment 1. 
In case of the ZZ — )■ e~''(l)e~(2)e'''(3)e~(4) process, for example, the function Fi, relative 
to the resonance assignment (1,2), (3,4) is chosen equal to 

Fi = {su [{si2 - Mlf + r|M|] X S34 [(S34 - Mlf + r|M|] }' , (A.l) 

while F2 is obtained by replacing 1^-7-3. In case of the ZW — )■ e^(l)e^(2)e^(3)z^e(4) 
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process, we have 



Fi 



{Sl2 [(Sl2 



Ml? + TlMl] 



[(S34 



Mlf + TlMl]]\ 



(A.2) 



for the region (1,2), (3, 4). For F2 one replaces again 1 -H- 3. In order to avoid possible 
errors in the analysis, the output at the level of Les Houches Interface is symmetrized again 
by randomly permuting the kinematics of the identical final state leptons. 

B. Born zeros 



All three production processes considered here have a Born cross section that vanishes in 
particular kinematic regions due to angular momentum conservation. In fact, since the 
vector interaction preserves chirality, the incoming quark-antiquark system has total spin 
equal to 1 in magnitude. If we consider the final state configuration where all leptons are 
parallel to the incoming quark, the whole system is symmetric under azimuthal rotations, 
and thus the angular momentum component along the collision axis equals the sum of 
the spins, and should be preserved. On the other hand, the final state leptons can be 
divided into two pairs with opposite helicity, thus yielding two systems with total spin 1 
(in absolute value) along the collision axis. Since the sum of the two spin 1 system yields 
either 2 or zero, angular momentum is violated in this configuration, and the amplitude 
must vanish. This property can be easily verified numerically. 

It was first pointed out in ref. [53] that problems can arise in the POWHEG formalism if 
the Born cross section vanishes or becomes particularly small in certain kinematic regions. 
We briefly recall the nature of this problem in the following. 

The basic POWHEG cross section formula can be written as 



da = B{^B)d^B 



A?„i 



^b) + ^1{^b) 



B{-^b) 



d^r 



+ [R{^) - R'{^)]d^, 



(B.l) 



where the real phase space is factorized in terms of the underlying Born ^b and the 
radiation ^^ phase space: d^ = d^rd^Bi and B, R and V are the Born, real and virtual 
amplitudes, i?* is an approximation to the real amplitude such that < R'^ < R, and that 
i?* — 7- i? in the limit of soft or collinear singularities. The choice R'^ = R is also allowed, 
and is often used. The variable t represents, in the present case, the transverse momentum 
of the radiated parton relative to the collision axis (thus t is a function of the phase space 
point), and to is a non-perturbative cutoff (of the order of a typical hadronic scale) on this 
variable. Furthermore 



and 



B{^B) = Bi<^B) + 



^l (^b) = exp 



Vi<^B) + / R'i<^)d<^r 



d^r 



R%^) 
B(<!>f 



(B.2) 



(B.3) 



lt>ti ^\^B) 

The Born term may become particularly small in certain kinematic regions; for example, if 
it vanishes for symmetry reasons. Besides the present case, another example of a behaviour 
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of this kind is the case of W production and decay, where there is a zero in the Born cross 
section if the outgoing lepton is anti-parahel to the incoming quark [53] , which is due to 
angular momentum conservation and to the left-handed nature of charged currents. It may 
then turn out that the B function does not vanish in the same region, due to the presence 
of the real term in eq. (B.2). From eq. (B.l), we see that, in this case, away from the 
Sudakov region, the real contribution may be enhanced by a factor B/B. The POWHEG 
BOX has a built-in mechanism to deal with this problem. If the flag withdamp is set, R^ 
is chosen to vanish in the regions where R differs too much (by more than a factor of 5) 
from its collinear or soft approximation, which are proportional to the underlying Born 
cross section. The contribution R — R'^, being non-singular, is then added independently. 
We thus conclude that the withdamp flag should also be activated in the processes we are 
considering. 

In practice, it turns out that the effect of the withdamp flag is much more important in 
the W'^W~ and in the WZ case, than in the ZZ case. This fact can be easily understood 
in the WZ case, where another region of vanishing Born cross section can be found which 
is less suppressed by phase space. Assume we have an incoming d quark in the positive 
rapidity direction, colliding with a u. In WZ production only the left handed component 
of the incoming d, and the right handed component of the incoming u, contribute, so that 
the incoming particles have a definite angular momentum projection along the collision 
axis equal to 1. Consider now the kinematic region where the W~ decay products are 
aligned along the collision axis, with the negative lepton in the negative rapidity direction. 
Again, the negative lepton is left handed, while the antineutrino is right handed. Thus the 
system has angular momentum projection along the collision axis equal to -1. Since the 
Z must also have zero transverse momentum, in order to balance the angular momentum 
component along the collision axis it should have spin 2, which is impossible. Thus this 
region is also suppressed, irrespective of the direction of the decay products of the Z. 

The W^W^ case is not as simple, since here there are two amplitudes that contribute: 
one with the two Ws emerging from the incoming fermion line (t-channel amplitude), 
and the other with the two Ws arising from the WWZ or WW^ vertex, the Z/^* being 
attached to the incoming quark line (s-channel amplitude). The same argument applied to 
the WZ case also applies to the t-channel amplitude. We can use a similar argument to 
show that the s-channel amplitude vanishes when the decay products of the two W^s travel 
along the same line, with the neutrino of one W is aligned to the lepton of the other. In 
this case, the decay system has angular momentum projection in the line of decay equal to 
2, and it can't therefore arise from a virtual Z or 7 decay. Although we cannot conclude 
that the full amplitude vanishes in either direction, even here we expect a strong kinematic 
suppression of the Born amplitude in certain regions. 

In Fig. 9, we show the transverse momentum distribution of the four lepton system in 
ZZ and W^W~ production if the withdamp flag is not set ^ (the effect in WZ production is 



^Because of the radiation zeroes, the upper bound of radiation tends to diverge as we increase the number 
of points used to compute it (i.e. the nubound variable in the powheg. input file). On the other hand, with 
a hmited number of points we still get a relatively small rate of upper bound violations in the generation 
of radiation, so that these distributions can be nevertheless computed. 
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Figure 9: Distribution oi ptAi for pp — J- W~^W~ — >■ e^v^iJ-^v^ and pp -^ ZZ -^ e'^e^fi^^', at 
centre of mass energy 7 TeV. The NLO results are shown in red, and LHE results using only the 
Born cross-section to generate events are shown in light blue. Also shown are the NLO LHE results 
with damping (in dark blue and labeled 'LHE-b') and without damping (in magenta and labeled 
'LHE-a'). The bottom plot shows the ratio of the NLO results to the LHE results with damping 
(dark blue, labeled 'b'), and the ratio of the NLO results to the LHE results without damping 
(magenta, labeled 'a'). 



very similar to that in W~^W~ production and is not shown). While we find a very modest 
effect in the ZZ case, in the W~^W~ case we find large enhancement of the differential 
cross section as a function of the LHE transverse momentum of the hardest parton. We 
also find that this anomalous behaviour disappears if we set B = B in PDWHEG (this is 
achieved by setting the flag bornonly to 1 in the powheg. input file). Thus, the anomalous 
growth is due to a large B/B ratio in some kinematical regions. 
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